Colloid-polymer mixtures in the protein limit 
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We computed the phase-separation behavior and effective interactions of colloid-polymer mixtures 
in the "protein limit" , where the polymer radius of gyration is much larger than the colloid radius. 
For ideal polymers, the critical colloidal packing fraction tends to zero, whereas for interacting 
polymers in a good solvent the behavior is governed by a universal binodal, implying a constant 
critical colloid packing fraction. In both systems the depletion interaction is not well described by 
effective pair potentials but requires the incorporation of many-body contributions. 

PACS numbers: 82.35.Np,61.25.Hq,82.70.Dd 



Adding polymers to suspensions of micro- and nano- 
particles induces depletion interactions that profoundly 
affect their physical properties. This phenomenon has 
important scientific and (bio-)technological applications. 
Polymers such as polyethylene glycol are routinely added 
to protein solutions to enable protein crystallization |l], 
U , a poorly understood process and of great importance 
in structural biology ||. In cell biology depletion inter- 
actions are key in the process of macro-molecular crowd- 
ing Q . Food and paint production are among the indus- 
trial sectors where depletion phenomena play a role. 

In this Letter we focus on mixtures of hard sphere 
(HS) colloids with a radius R c and non-adsorbing poly- 
mers with a radius of gyration R g , in the regime where 
q = Rg/R c > 1. This is often called the nano-particle or 
"protein limit", because in practice small particles such 
as proteins or micelles are needed to achieve large size- 
ratios q. Whereas the opposite "colloid limit" (q < I) 
has been well studied, the physics in the protein limit is 
less established. This imbalance is partially due to the 
lack of well characterized experimental model systems 
for the protein limit and partially to a poor theoretical 
understanding. The colloid limit can be well described 
within the framework of effective depletion pair poten- 
tials H q|, in contrast to the protein limit, where the 
interactions cannot be reduced to a pairwise form |t], [|. 
Nevertheless, for biological and industrial applications, 
this regime is at least as important as the colloid limit. 

One of the first theoretical treatments of colloid- 
polymer mixtures in the protein limit was by de 
Gennes ||, who showed that the insertion free-energy 
F c (1) of a single hard, non-adsorbing sphere into an ather- 
mal polymer solution scales as 
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when R c < 
€(<f>p) ~ R gVp 

1/fcgT is the reciprocal temperature, 
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Here, /3 — 
0.59 is the 



Flory exponent and 



Pp^knRq is the polymer vol- 



ume fraction for a polymer number density p p , so that 



(f>p « 1 at the crossover from a dilute to a semi-dilute solu- 
tion |l(J. The prefactors can be calculated by the renor- 
malization group (RG) theory |Q, yielding: /3F C (1) w 
4.390 p (7 -1 ' 3 which has been verified by computer sim- 
ulations for small q [p^f . Based on this description of 

Fc , de Gennes |l(| and Odijk || predicted extensive 
miscibility for colloid-polymer mixtures in the large q 
limit if R c < £. However, it is well known that protein- 
polymer mixtures do phase-separate (l3[ Q. Recently, 
Odijk et al. |l5| suggested that a poor solvent could fa- 
cilitate phase-separation. Sear |l(| altered the form of 

Fc to include effects when R c » £, and also predicted 
phase-separation with a truncated virial theory. The 
same author recently proposed an alternative theory |l7j 
where the colloids induce depletion attractions between 
the polymers, leading to a poorer effective solvent and 
eventually phase-separation. Mean-field cell model cal- 
culations also predict demixing (l8| . Another promising 
approach uses integral equation techniques fl9|| to predict 
spinodal curves and critical points. However, all these 
theories suffer from several uncontrolled approximations 
leading to different predictions for the causes and prop- 
erties of the phase-separation. To clarify this situation, 
we performed computer simulations with as few simpli- 
fying assumptions as possible, on which we report in this 
Letter. 

We have recently used a coarse-graining technique p3] 
to study the colloid limit, and found quantitative agree- 
ment with experimental fluid-fluid binodals pl| , and sig- 
nificant qualitative differences between interacting (IP) 
and non-interacting (NIP) polymers. Here, we study the 
same athermal model of HS colloids and non-adsorbing 
polymers in the protein limit, and calculate, for the first 
time, the full fluid-fluid binodals by direct simulation. 
The results for the IP and NIP show even larger quali- 
tative differences, and many-body depletion interactions 
must be invoked to understand the phase-behavior. 

The simulation model consists of polymers on a sim- 
ple cubic lattice mixed with HS colloids. The interacting 
polymers in a good solvent are modeled as self avoid- 
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FIG. 1: Fluid-fluid binodals for a mixture of non-interacting 
polymers and HS colloids with different size-ratios q. Crosses 
indicate the estimated critical point, obtained by extrapolat- 
ing the calculated phase boundaries. The full lines are a guide 
to the eye. Dashed lines denote the simple theory described 
in the text, with stars marking the critical points. 



ing walks (SAW) of length L, which have a radius of 
gyration R g ~ L u . The non-interacting polymers are 
modeled as random- walks, for which R g ~ L 5 . In both 
models there is an excluded-volume interaction between 
the colloidal HS and the polymer segments. The sim- 
ulations were performed on a D 3 lattice with periodic 
boundary conditions, where D — 48 and D = 100 for 
the NIP and IP system, respectively. Throughout this 
Letter, we use the lattice spacing as the unit of length. 
For the NIP the colloidal HS diameter was a c = 5.5 and 
the polymer length was L — 50, 100, 200 and 500, cor- 
responding to q = 1.03,1.45,2.05 and 3.2, respectively. 
For the IP L = 2000, and cr c = 10, 14 and 20, yielding 
q = 3.86, 5.58 and 7.78, respectively. Colloidal positions 
have continuous values, but when we calculated the in- 
teraction between colloid and polymer the colloids were 
shifted such that they occupied a constant number of 
lattice sites to prevent spurious attractive positions for 
single colloids (other lattice effects, although unavoid- 
able, are expected to be small.) Thermodynamic state 
points were calculated in the grand-canonical ensemble, 
i.e. at fixed volume V, colloid chemical potential fi c and 
polymer chemical potential /i p using Monte Carlo (MC) 
techniques. The NIP were sampled using an (exact) lat- 
tice propagation method |^2|, ^3) , while the IP configura- 
tions were generated using translation, pivot moves and 
configurational bias MC pi| in an expanded ensemble to 
facilitate insertion of long chains p5[ . Typical simula- 
tions lengths are 10 9 Monte Carlo moves per state point. 
In order to determine the liquid-liquid binodals we first 
estimated the coexistence line by scanning a series of p, c 
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FIG. 2: Fluid-fluid binodals for a mixture of interacting poly- 
mers and HS colloids at different size-ratios q. Filled symbols 
are direct simulation data. The open symbols are the colloid 
limit (q < 1) results from Ref. [Ell. Solid lines are a guide to 
the eye. Inset: The same binodals plotted in a reduced poly- 
mer density representation. The dotted curve corresponds to 
a simple theory for the universal binodal when polymers are 
in a good solvent while the dashed line is for polymers in a 
poorer solvent. Crosses and stars as in Fig. m. 



for several values of fi p and locate the fi c for which a sud- 
den density change occurred. Subsequently 8-10 (fi p , /i c ) 
coexistence state points were simulated simultaneously 
using parallel tempering [p5[ . When the estimated co- 
existence points are sufficiently close to the true binodal 
and to each other, and near the critical point, this scheme 
results in proper ergodic sampling of both phases. If 
necessary, the chemical potentials were adjusted towards 
coexistence. We used the multiple histogram reweight- 
ing ^(| technique to determine the precise location 
of the (fi c , n P ) coexistence line, and the phase boundaries 
in the (</> c , </> p ) plane, where <p c = p c ^irR 3 is the colloid 
volume fraction, with p c the colloid number density. 

Figs. [j] and || contrast the calculated phase diagrams 
for NIP and IP for several size-ratios q. Firstly, we note 
that both models show extensive immiscibility, in agree- 
ment with experiment Jf4"| . Secondly, the two systems ex- 
hibit striking differences: for the NIP, the critical colloid 
volume fraction c/>£ tends to zero with increasing size ra- 
tio q, while the IP exhibit a nearly constant value of 
For both systems the critical polymer concentration <\> c v r%t 
increases with increasing q. The phase-separation occurs 
well into the semi-dilute regime for the IP, again in qual- 
itative agreement with experiment Jl3| . Properties of 
semi-dilute polymer solutions are independent of poly- 
mer length, being instead determined by the correlation 
length £, which is a function of the monomer density 
c = Lp p . The phase behavior of the polymer-HS mixture 
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should therefore only be a function of the ratio i? c /£ [ P^I - 
Indeed, when the phase lines in Fig. |2] are rescaled with an 
accurate prescription for £(p p ) Q, the binodals nearly 
collapse onto a "universal binodal" , as shown in the inset 
of Fig. |^. This explains why the critical colloid packing 
fraction is nearly constant in the simulations. Similarly, 
4> p rlt scales as ~ g( 3l/ ~ 1 )/ iy pa q 1 - 3 . For comparison, 
we have also included results for q < 1 from Ref. |2l[] 
in Fig. pi These results do not exhibit the same scaling 
behavior, since they are not in the semi-dilute regime. 

The differences between NIP and IP phase behavior 
can be rationalized with some simple theories. Consider 
a Helmholtz free-energy F of the form f3F/V = f = 
Ic S + fp + fcp- Here, the HS free energy f^ s is given 
by the accurate Carnahan-Starling expression J27|, and 
the polymer free energy L, for either IP or NIP solu- 
tions is well understood p0[ . The contribution due to 
the HS-polymer interactions f cp is non-trivial. A first 
approximation truncates after the second cross-virial co- 
efficient, yielding f cp ~ p c Fi X \ For NIP the insertion 
free-energy F^ is exactly known so that f cp takes 
the form f%* = p p (f> c (l + ^| + 3g 2 ) = p p (j) c b cp , which 

defines the reduced cross-virial coefficient b cp . Since fl p 
grows with increasing q, immiscibility sets in at lower col- 
loid packing fraction <p c . The theory can be improved by 
realizing that the polymers only exist in the free volume 
left by the colloids j28|. Simply taking this free volume 
to be 1 — 4> c is an adequate first approximation for the 
protein limit. The trends for the binodal lines calculated 
from this simple theory, shown in Fig. [l], agree quali- 
tatively with the simulations. For example, the critical 
point shifts to smaller (f> c and larger <j) p for increasing q, 
and the binodal lines cross at a low cj) c . For computa- 
tional reasons the simulations only go up to q = 3.2 and 
we expect better quantitative agreement for larger q since 
^crit d ecreases so that the second-virial theory should 
become more accurate. In the q — > oo limit, this theory 
yields <f>f u -> l/"b cp ~ l/(3(? 2 ), and <f$* = q 3 /b cp ~ q/'S. 
Note that in the same limit, the penetrable sphere or 
Asakura-Oosawa model |28| scales somewhat differently: 
^crit ^ i/q3 anc j ^crit ^ g ear g a i re ady pointed out 

the <\fj" it — y behavior using a slightly different prescrip- 
tion for the free volume than we employ here. Here we 
claim that the limiting results are a general feature of 
free- volume theories. 

In the IP case, f cp is more difficult to estimate, even 
for a second cross-virial theory. The R c <C £ limit is 
given by Eq. (|l|) with the prefactors from RG theory. 
For R c > £ we have previously shown that Fc is given 
by F c (1) = §7ri?3n + 6irR 2 cls @, where the polymer os- 
motic pressure n ~ is well known p0| . However, 
since Eq. (|f|) is essentially a surface (depletion layer) 
contribution, we use a simple approximate second cross- 
virial term f cp = p c ((3H(p) ^tt R 3 , + 4.390 p q -1 - 3 ) , which 



reduces to the correct form in both the R c <C ^ and the 
R c S> i limit [^3| . As with our treatment of NIP, we take 
the effect of the colloid excluded volume into account by 
computing the polymer densities in the free volume frac- 
tion 1 — (f> c (see Ref. |2^] for a complimentary approach). 
The theoretical binodals were calculated using accurate 
expressions for £ and n JI2[ and are compared with the 
simulation results in Fig. in the R c /£, versus <fi c plane. 
The qualitative agreement suggests that we can also use 
this theory to estimate the effect of a poorer solvent on 
the binodals. Following Ref. Jl5| ], we alter the scaling of 

£ t° f ~ 4>p <5//3 so th at n ~ ~ p p 0p _1 . Interestingly, 
Fig. H shows that using S ~ 1.5 instead of the appropriate 
exponent for polymers in a good solvent (<5 ~ 2.3), does 
not result in important differences in the binodals. Of 
course, the differences will appear larger in the (</> c , (j) p ) 
plane due to the different scaling of £. One must keep 
in mind, however, that these predictions follow from a 
simple scaling theory and qualitatively different behav- 
ior may emerge when one approaches the 0-point (where 
5 = 1). 

To illustrate the many-body nature of the depletion 
interaction we estimated the phase behavior by approx- 
imating the system by colloids interacting via pairwise 
effective potentials. We computed the effective pair in- 
teraction v(r) between two colloids in a bath of IP's, by 
j3v(r) = — \ng(r) for p c — > 0. The colloid radial dis- 
tribution function g(r) was estimated by measuring the 
insertion probability of a HS at a distance r from a sec- 
ond fixed HS in a SAW polymer solution using the above 
MC techniques. Results for a single size-ratio q = 7.78 as 
a function of (f> p are shown in Fig. [| Several features are 




FIG. 3: Effective colloid-colloid pair potentials induced 
by interacting polymers for q=7.78. Theoretical lines from 
Ref. pc(]. Inset: Reduced second osmotic virial coefficient 
as a function polymer densities for several 

size ratio's. 
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similar to the colloid limit j30|: the range shortens and 
the well-depth increases with increasing <p p . Interestingly, 
our simple depletion potential Q, derived for the col- 
loidal limit, also works semi-quantitatively in this regime. 
A good measure for the attractive strength of an effec- 
tive pair potentials is given by the second osmotic virial 
coefficient p7| , shown in Fig. ||. The saturation of £>J 
for larger q is an interesting qualitative feature: appar- 
ently the shortening of the range compensates the deep- 
ening of the attraction, so that the total cohesion does 
not increase with increasing <p p , something also found in 
RG (ll]] and integral equation calculations 1 19 . For pair- 
wise interacting systems, phase-separation typically sets 
in when B^ < —1.5 [[n). Here, the saturation of sug- 
gests that for large q the pair interactions do not provide 
enough cohesion to explain the phase-separation. We ar- 
rive at the same conclusions with simple mean-field theo- 
ries H?! , which should be relatively trustworthy given the 
long range of the pair potentials. Obviously, for g > 1 
a pair- level description is not sufficient, and many-body 
interactions must be invoked. 

For the NIP, one might also expect many-body inter- 
actions to be important for large q. A good approxima- 
tion to the pair-potentials exists |2^, ||(J , from which the 
second virial coefficients at the calculated critical points 
follow: B%(q = 1.03) w -13.1; B$(q = 1.45) « -16.4; 
B%(q = 2.05) w -22.7. Even though the actual crit- 
ical </> c 's are very low, so that a second- virial descrip- 
tion might be thought to be sufficient, the analysis 
above shows that for NIP the pair-interactions provide 
too much cohesion, opposite to the IP case. Clearly, 
many-body interactions must also be invoked to describe 
the phase-behavior correctly, as suggested by other au- 
thors gig §3, 



and R. Tuinier for helpful discussions. 



In conclusion, we have shown by computer simulations 
that a mixture of polymers and non-adsorbing HS col- 
loids shows extensive immiscibility in the protein limit, 
where the polymer-colloid size-ratio g > 1. For IP the 
phase-behavior is dictated by a universal binodal in the 
semi-dilute regime. For NIP, the colloid packing frac- 
tion tends to zero for increasing polymer length. In con- 
trast to the better studied colloid limit, pair interactions 
are not sufficient to rationalize the phase behavior. We 
hope that future experiments on HS colloids with non- 
adsorbing polymer will test these predictions. Future 
work might include extensions to non-spherical particles, 
poor solvents, and adsorbing systems. 
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